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^S) ■ In the present paper an approach for investigation of the dis- 

t*^ ' ordered two-component Ising systems with long range interaction 

C""~~ , has been suggested. Possible applications to metalic and magnetic 

0^ ■ alloys and lattice gas are considered. We have also obtained nu- 

merical results for thermodynamical properties of these models. 
The comparison of numerical results obtained within mean field, 
gaussian field and two-tail approximations are carried out. 
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1 Introduction 

The present paper is devoted to the investigation of a binary annealed 
disordered magnet on the basis of Ising model. This model and its sim- 
pler versions are often used in the study of metalic alloys H, g, ^ , mag- 
netic alloys [Q g|, ferroelectric solutions, in particular, the ferroelectrics 
with hydrogen bonds [|6| . 

Let us point out several works concerning the magnetic alloys. In 
the work jj], V.G.Vaks and N.E.Zcin obtained phase diagrams for the 
model of the binary magnetic alloy with non- magnetic impurities within 
two particle cluster approximation at T > Tm {Tm is the temperature of 
the magnetic ordering). It was noted that in mean field approximation 
(MFA) in contrast to cluster approximation the magnetic subsystem does 
not influence the phase diagram. Later by means of MFA T.Kawasaki g 
investigated the influence of the magnetic subsystem on the properties 
of the atomic subsystem and v. v. at T < T,„. In the work B a detailed 
investigation of this model within two-site cluster approximation were 
performed. There was found the influence of the spin subsystem on bin- 
odal (coexistence temperature Tb) and spinodal (spontaneous separation 
temperature Ts) temperatures of the system in a spin ordering phase. 
Here the difference between quenched and annealed types of disorder has 
been investigated and correlation functions have been calculated. 

The non-magnetic version of the model is often applied to the study of 
real binary alloys El 0, H] . Using the pseudopotential theory for defining 
of intersite interaction Z.Gurskii and Yu.Khokhlov tried to explain the 
properties of Ca-Ba and K-Cs alloys B] . For the thermodynamical aver- 
aging they used the collective variables method with a non-selfconsistent 
equation for unknown chemical potential. 

This work has been directed to the study of a disordered model with 
long range interaction or short range one in a case of large z [z is the 
number of nearest neighbours). At first we consider the Hamiltonian 
of A/-component Ising model and study the model within mean field 
approximation (MFA), Gaussian field approximation (GFA) and two-tail 
approximation (TTA). But for the sake of simplicity all approximations 
are formulated for the case M = 2. 



2 Description of the model 

Let us consider M-componcnt system with site disorder on Bravais lat- 
tice with the foUowing Haniiltonian: 

-(3H = H = n^ + ns, (2.1) 



(2.2) 



where 




H, 


M N M N 
a=li=l a,l3=lij = l 


Hs 


M N M N 
a=li=l a,l3=li,j = l 



All parameters of the Hamiltonian (2.1) contain the factor (— /3). The 
Hamiltonian Ti^ corresponds to the disordered ionic subsystem. Here 
£q = ^q + /^a, /^q = ~Pl^a i'' the dimcnslonless chemical potential of 
a-type ions, "da is related to the difference between alloy component 
ion characteristics; Vafj{i — j) = Va/siRij) is the potential of effective 
ionic interaction; Xia = 1, if the site is occupied by the ion of type a 
and otherwise. The Hamiltonian TCs corresponds to the system of N 
Ising spins {Sia — Sf^) situated on the ions with external field energy 
ha — ifJ'oM) and exchange interaction Tap(i — j) = lap^Rij)- 

In the present paper we confine ourselves by the case of binary al- 
loy (M = 2). It should be noted that suggested approximations below 
can be applied to magnetic alloys with arbitrary number of sorts. It is 
convenient to use the spin variables instead of x^i, 2:^2(0: ~ 1, 2) 

_ 1 ^ (^l)"'5'iO CI _ (0 A\ 

Xia — n J ^iO —■ Xii — Xi2. V^-^j 

Now the Hamiltonian (2.1) can be written as follows 

n^s = E{e)+nQ{T), (2.5) 

N -. N 

UeiT) = X Xr-0«.- + 2 51 X I^.<r'{R^J)Q^.QJa■ 
(J— 0,1,2 i— 1 a. a' i.j — 1 

Here we use the following notations 





I{R) 

I{R)= I Xii(i?) Ji2(i?) 

Xl2{R) l22{R) 

j^E = i(ei + £2) + ^(^11 + V^22 + 2T/12), (2.6) 

To - ^(£1 - S2) + i(Vll - F22) (y„,/3 = K.,/3(g ^ 0)), 

4J(i?) = T/ii(i?) + V22iR) - 2Vi2{R). 

The Hamiltonian 7-^e(r) describes the three-sort spin system. The ther- 
modynamical potential (TP) of system (2.1) has the foUowing form (h = 
{hi,h2)) 

f^(£,h)-~/3-iln%,,,}e^(^'^)=S(e) + f7e(r). (2.7) 

Here r2e(r) is the TP for the Hamiltonian Hq^T). The correlation 
functions (CFs) of an arbitrary order can be obtained on the basis of 
the relations 

n^i\z^s,, ..., ii5i) - (e,,5, • • • e„5,)' = j^ -^[-/3f^e(r)]. 

(2.8) 
Here we have replaced the uniform fields Fq with nonuniform ones Via- 
The TP (2.7) is a function of temperature (T) and chemical potentials 
(/Ui,/Z2). Transition to concentration variables Ca instead of fia can be 
performed by the Legendre transformation 

^F{c, h) = lo(£, h) - ^ ^iy^. (2.9) 

Q=l,2 

Here -F(c, h) is a free energy of system (2.1) and chemical potentials are 
defined by (too = {Sm)) 

'-t^£i-'-(..,i^i^i^^ (.10, 

The system (2.10) gives only one independent equation for /ii — ^2 or Pq, 
because ^^ Ca = 1- Here we consider the system which undergoes both 
the separation phase transition and magnetic phase at temperature Tm- 
Therefore, the expression for {Sio) can contain the magnetic moments 
Til — (Sii) and m2 = ('5'i2) as well. The equations for them we will 
obtain proceeding from (2.8): 

mi = ^ • ^hpneiT)], "^2 = ^ ■ ^Mf7e(r)]. (2.11) 



The system of equations (2.10), (2.11) gives the values for To, mi, m2 as 
function of temperature T and concentration c = ci . The spontaneous 
separation transition occurs at Tg (spinodal temperature) which is found 
as divergence temperature of CF ((S'io<S'jo)^)q=o- The coexistence tem- 
perature Tb (binodal temperature) of both phases A and B is calculated 
from the following system 

rj(e(mo), h, Tb)/mo=m„A = ^^(^(too), h, Tb)/mo=m„B (2 12) 

To{moA,h,Tb) ^Ta{moB,i^,Tb) 

The system (2.12) gives the dependences moA{Tb),moB{Tb) and it is 
solved jointly with system (2.10), (2.11). 

Although for the system (2.12) exact solution has not been found, 
one can make some conclusions about the symmetry of phase diagrams. 
In case Xn = X22 and Fi = F2 the thermodynamical potential 



f^e(ro,h) = -/3"4nSpexp 



N 



1 ^ 



y^Fo(xa-a:i2)+ 

.r!e(-Fo,h) 



will be even function of Fq (because trace over xn and Xi2 becomes 
equivalent). Similarly, it can be proved that 

TOo(ro) = {Xil-Xi2) = 
Sn(xii-a;i2)e5I)'"''*^^"~^*2^+2Z]^(^^-')(^'i^^*2)(2:jl-2:i2)+'Hs 

SDeS'"'''-^"^^'^-'^^'^'-*^-'-"-^"^^'^-''^^^^^^^-'^^^ 

Sp(a;-2— Xi)eS''^'""-'''^"^'^'^''^^5I]-^(*^-')('^»i^^*2)(2;3i-a;j2)+'Hs 

= {xi2-Xii) == -mo(-ro) 

is odd function of Fg, where Fq ^ {pi — ^2) and determined by the 
formula (2.6). This means that f^e is even function relatively variable 
Too and first equation of the system (2.12) has following solution mos = 
—rriQA- Then second equation for ttiq can be rewritten as following: 
ro(^^Oyij h, Tb) = 0. Therefore, phase diagram is symmetrical relatively 
variables ttiq or permutation of components' concentration ci and C2. 

The spinodal temperature also is symmetrical function of ttio = ci — C2 
in case Tn — X22 and Fi = F2. Spinodal decay takes place if the CF 



{{SioSjo)'^)q=o diverges. It means that the system becomes unstable 
against the infinitesimal fluctuation of concentration: -f-?- = 0. Since 
■S^ = Tn , the last condition may be rewritten as -S^ = 0. Fn is odd 
function of mo, then ^—^ will be even function of rriQ. This means if 
for a certain temperature Tg the equation 4-^ — has solution toq , the 
solution — TTiQ also exists and spinodal temperature Ts is an even function 
relatively rriQ. 

Further for calculation of TP and CFs of system with Hamiltonian 
T-Cq{T) we will use the expansion over 1/z and the results of ref. 0. 
Now the TP J^e(r) can be written as follows {Tafi = (Iq/sIq ~^ Q),V — 

-|:ne(r) = -i Y. ^ss'msms, + fis,), (2.13) 

<5,<5'=0,1,2 

/(ae) = (F(")(ae + ^))-i ^ \f^, {F^^} (33 + a)) ^ 

S,5'=Q,1,2 

-^El^det[l-i(<z)(Fi,',^(ae + ^))]. (2.14) 

Here irreducible part /(ae) of r2e(F) is written with the accuracy up 
one-loop diagrams (two-tail diagrams approximation (TTA)). 

We use the following notations from ref. |Q for averages of arbitrary 
function y{se + cr) over fluctuating fields cr with Gaussian distribution 
function /92(c): 

(z/(ae + cr))= daod<7ida2P2{(T)y{se + cr), (2.15) 

p,{a) = [dct2^A(2)]-i/2cxp{-i^[A(2)]->,a,,}, 

^ 5,5' 

aeo = Fo +Imo, Xa == h^ + Ylaf3mp {a, P = 1,2), 

/3 

where [A*^^^]"^ is the inverse matrix to A^^^ . In (2.15) we use the following 
notations for functions on x = [xq, xi, X2) 

F(°)(x) = InZ(x) = ln{e"''Zi(a;i) + e-"°Z2(a;2)}; 

^i"(x, = A^.«.(x);^£),x,.J--A^<«.(x);P.16) 

Za{xa) = Sps^e^"^" =^ 2cosha;o,; 

Sq=±1 



From stationarity conditions with respect to nine variables rriQ^ma, 
-^00 ' -^Oa ~ '^qO' ■^a0 ^^ '^^'^ ^^^ system of nine equations for nine 
unknown Tcto^, a{,o\ A^,^^ = A^^J, A^^j 

2c-l = mo-(fo^'^(ae + ,T)}, m„ = (F«(ae + ,t)}, (2.17) 

xfl = [^Y.'^i-imF^'\^+cT))]-^im . (2.18) 



In the works of Onyszkiewicz at aU ||] GFA was suggested. It can be ob- 
tained expanding lndet[i-X(g)(FJ5^?(ae + cr))] in (2.14) and i-^ EJl- 

J(g)(/'(2)(ae + cr))]-ii(q)) in (2.18) in 2(17) up to the terms of second 
order. In the result : 



/(ae) = (f(")(ae + <T))-i Y. 4''(^i'' (^ + '^)) ^ 

<5,5'=0,1,2 

-^ E (^ifi'(^ + '^))(<1(«' + '^)) E^^^^(9)I.'..(9),(2.19) 

i5,(5i, (52=0, 1,2 g 

^i'-' - ^ E (^i,t(ae + <T))5]j,,,(q)J,,,,(g).(2.20) 

(5i, (52=04, 2 <? 

In contrast to GFA and TTA MFA permits to get analytical re- 
sults for thermodynamical functions and transition temperatures in case 
of non-magnetic alloys. This approximation one can obtain neglecting 
fluctuation of molecular field rtif, . The expression for TP is 

--^r!e(r) = -i Y. l5S'msms'+F<'''\as). (2.21) 

<5,(5'=Oa,2 

The equation for chemical potential and order parameters is determined 
by the relations (2.10), (2.11). In MFA one can obtain: 

2ci-l=TOo=^o^''(ae)= [e="°^i(a5i)-c-==^"Z2(ac2)] ^(ae)"'(2.22) 
m„ - F<(i)(ae) - e(-i)°^'-°Z«(ae„)Z(ae)-\ (2.23) 

The equation for ttiq can be transformed to two dependent equations for 
the components' concentration Ca (see 2.10)): 

c„ = e(-l)°^''^°Z„(ae„)^(ae)"^ (2.24) 



From (2. 23), (2. 24) it follows 

m„ = c„^(i)(ae„). (2.25) 

The equation (2.22) for unknown difference of chemical potential can 
be solved analytically 

To^^Imo + -\n- -— — -. 2.26 

2 1 — mo Zi(Xi) 

Using (2.9) one can obtain following result for the free energy: 

^F(c,h)= J2 Ca^a+ln2 + ^(T4i+ 1/22 + 2^12) 

+ Ta^22 - ^11 mo + -Iml - - In 1 - m^ - -^ In ^ 

4 2 2 2 1 — niQ 

-]- V I„^^m„m^-^ln|4^ + ilnZi(aei)Z2(ae2). (2.27) 

The system (2.12) in a case In = 222 and Fi = F2 can be solved 
analytically with the following result: 

T, = -'^7^ ^ . (2.28) 

1+mo Zi(aoi) 

However, agi and ae2 is to find from equations (2.11). The expression for 
spinodal temperature can be found exactly only for non-magnetic case: 

Ts =I{1- ml) ^^ciC2. (2.29) 

The pair CFs we obtain on the basis of the relation (2.8) (the de- 
tails see in [Q) and this procedure leads to expression for matrix of CF 



((a^^ia^ji)") -^^e,oo(q)' (i^^i^jcSjo:)'') =^^ffl^{ci^-iQ) 



Here /] l^ (q) is the irreducible part of CF (Iq (q) with respect to in- 
teraction X(q). It is obtained by double differentiation of irreducible 



part/(ae) with respect to nonuniform fields sSi^Si, iBi^s^- Within TTA at 
T > Tm Z^^-* is independent of q and has form 

/i'i(q) = (Hli^ + -)> "^ ^i'l(^), (2.31) 

p(2). ^ 2i''(aOc)e \MFA (2.32) 

^ (2) 

From matrix expression for ^q (q) (2.30) it can be found 

<oo(q) = [1 - $oo(q)T^(q)]-'*oo(q), 

f}[|)^^(q) = [1 - <l^(q)J^(q)]-l<l^(q), (2.33) 



$oo(q)=,/,io'(q)+ E .C(q)(i^(q)[i-r(q)i^(q)]-M /^'l,(q), 

aia2 

Kc'ici) = /l'i(q) + /lo'(q)^(q)[l - /oo(q)(q)]-\ff,(q), 

r(q) = {/S(q)}' I^(q) - {laMq)}- (2.34) 

At the temperature of second order phase transition the CFs of system 
diverge at certain values q* (we consider here only case q* = 0). 

In the case T > T^ the /^^\q) = and $oo(q) = /(i?(q), *^„'(q) = 

(2) 

/aa'(q)- Within MFA from first expression we obtain the Tg (the in- 
fluence of magnetic subsystem is absent) and from the second one the 
equation for Tm {Tm is the temperature of magnetic transition) 

Ts = 4ciC2Z(0), 

2Tm - ciTn + CT22 - {(ciTn - 01^22)' + 4ciC2 T^2V/\2.35) 

To/j^Ia/s/^F^. (2.36) 

3 Results of numerical investigation 

The Hamiltonian (2.1) corresponds to the most generalized model of 
M-component magnet. There is no possibility at present to obtain any 
analytical solutions in this case. Therefore, in this section we represent 
some numerical results performed for particular cases of one-component 
Ising model, binary alloys and lattice gas. Here the numerical calcula- 
tions become simpler. 



3.1 Ising model 

The case of Ising model will correspond to M = 1. The ionic part of the 
Hamiltonian becomes constant and does not influence on the thermody- 
namical and correlation functions of the system: 



N ^ N 

7i = 2J Ti^ii + X 2J ^uSiiSji + const (3.1) 

i—l ^,i — 1 

For calculation of free energy, magnetization, two-tail we can use 
formulae (2.13)-(2.20) and put a, f3 ~ 1. In this case we have to solve 

system of only two equations for order parameter mi and fluctuating 

(2) 
parameter Xn ■ Therefore, one can write the following expression for 

generating function in the framework of two-tail approximation (2.13)- 

(2.18) that corresponds to free energy in the case of Ising model: 



-^EHl-^ii('Z)(i^i'i'(aci+'Ti))]. (3.2) 

The equations for order parameters mi and two-tail \\-( follow from the 
equations (2.17)-(2.18) : 

mi^{F[^\sii+ai}), (3.3) 

Xf) = ^y hl^ . (3.4) 

Here the averaging is performed with distribution function : 

P2{a} = /f e<-e-K^Af) ^ ^_,~^ _ (3.5) 

GFA gives the following results for the free energy and equation for two- 
tail \f} (see (2.19), (2.20)): 

~J^F{V) = -\liiml + (F(0)(8ei + ai)) Ia^^ (F^f (aci + ^1)) 



+ \{F[f{^i+ai)) 



(^i^?^(«=i+'^i))^E^ii(«)-24i^ 



9 



(3.6) 



mi = (F«{aei+ai})^^, (3.7) 

Afi^-(^i'i^{«^i+^i})p,^E^ii('^)- (3.8) 

q 

Neglecting the fluctuation term —^^Iij/S.Si/S.Sj in the Hamiltonian 

ij 

(2) 

we get the MFA approximation. This concerns the case \\i = 0, 
-(3F{V)/N = -\liiml + F(o){aei} and mi = F^^^{sii]. 

The numerical calculations were performed for the model with near- 
est neighbours interaction : 

3 

Xii(g) = /yj cosga. (3.9) 

a=l 

The results of this investigation can be seen in Fig.l, 2. In Fig.l the 
free energy of the system as a function of order parameter for various 
temperatures is depicted. It is known, that we get first order phase 
transition in the framework of two-tail approximation. This can be seen 
in Fig.l. The curve for T/I{0) ~ 0.72 has two minima with equal energy. 
Although the GFA is more rough than two-tail approximation, it gives 
us a second order phase transition at T/I{Q) — 0.86. The results of 
temperature dependences of magnetization in different approximation 
are plotted in Fig. 2. 

The obtained results show that the two-tail approximation predicts a 
non-physical behaviour in the close vicinity of the critical temperature. 
This may lead to the non-controlled errors by investigation of binary 
alloys and lattice gas. Therefore, for these systems we shall restrict only 
GFA and compare results with MFA. 

3.2 Binary alloy 

The Hamiltonian of the non-magnetic binary alloy follows from (2.1) if 
we put Iii(i - j) = Ji2(i - i) = l22{i - j) = 0, Fii = ri2 = 0: 

2 N 2 N 

a—1 i—1 a,/3— 1 i,j — l 

The MFA and GFA may be easily reformulated for this case using formu- 
lae of previous section. Our numerical calculations concern the system 



10 



-0.53- 




-0.57- 



-0.61- 



-0.65 



T = 0.86 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1: The free energy -^F as a. function of magnetization m for 
various temperature (solid curves - two-tail approximation, long dashed 
curves - GFA, short dashed curves - MFA) 
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Figure 2: The magnetization m as a function of temperature T (solid 
curve - two-tail approximation, long dashed curve - GFA, short dashed 
curve - MFA) 
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Figure 3: The two tail A as a function of temperature T (solid curve 
two-tail approximation, dashed curve - GFA) 
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with nearest neighbours interaction. As may be proved in section 2 
the critical temperatures of binodal and spinodal decay depends only on 
1- ~ Vii + T^22 — 2Vi2 and are symetrical relatively substitution C2 instead 
of c\. This can be seen in Fig. 4, where binodal and spinodal temper- 
atures as a function of c\ are depicted. The comparison of GFA and 
MFA approximation results show that GFA approximation is essentially 
better that MFA in the vicinity of C\ = 0.5. When ci ^ or ci — > 1, the 
difference disappears. It should be noticed that two-tail approximation 
will give non-physical results namely in the vicinity of c\ = 0.5. 

3.3 Lattice gas 

We shall get the Hamiltonian of the non- magnetic lattice gas, if we put 
^12 (^ — j) = V22(i — j) = in case of binary alloy: 

N N 

']-l^^^Xi + -^V{i-j)xiXj. (3.11) 

2 — 1 ^J — 1 

Indeed, the Hamiltonian (3.11) may describe the gas. Let us consider the 
system in volume V and divide volume on A^ = V equal cubic cells. We 
will use the condition that more than one particle can not be situated in 
one site. Therefore, we get a lough model of a gas with repulsion between 
particles on short distance and model interaction for long distance. 

Here we neglect the kinetic energy of the particles. But being taken 
into account it leads only to the renormalization of the chemical potential 

^ ^ ^ + - In (27r/i" VfcT) , (3.12) 

where m - particle's mass. 

Similarly to (2.5), the Hamiltonian (3.13) can be transformed to 

N N 

n = E{p) + Y, To^^o + 2 H ^(* - ])SmS,o, (3.13) 

i—1 ij — 'i- 

where Eip) = i^ + i^, To = ^^i + \V , 4J(* - j) = V{i - j). 
It is known from thermodynamics, that thermodynamical potential of 
the system fl{^, T) = —pV, where p is pressure and V is volume of the 
system. The average of Xi that equal n = 1/v, plays the role of gas 
density. The state equation can be written in the following form 

pV = r!(/i(n), T) orp^ ^n{^i{n), T). (3.14) 
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Figure 4: Phase diagrams for the model (3.9): 1- correspond to the 
binodal temperature, 2- spinodal temperature (sohd Une - GFA, dashed 
hue - MFA) 
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In MFA the state equation can be obtained in the analytical form. 
The expression for TP follows from (2.21) and has the following form: 

^niii,T) = ^^i + ^V-^Jml + ln2cosh{To+Imo),i3.15) 

where ro(mo) = -/Jlmo + ^In^^, ^ = 2ro - ^I. If we take into 
account that mo = 2n — 1, we can obtain state equation in variables p, 
n, T: 

- = -/3an2 + $(n), (3.16) 

where $(n) = In-r^- It should be noted that for van der Waals state 
equation <I>(n) = j^j^- 

We have also obtained isoterms in GFA. Since the Hamiltonian of 
the lattice gas is identical to the Hamiltonian of the binary alloy with 
Vi2ii — j) = V22{i — j) = 0, we can use all approximation formulated in 
the previous subsection to its investigation. 

In the Fig. 5 the isoterms for the lattice gas with nearest neighbours 
interaction are depicted. The dashed curve corresponds to the MFA and 
solid curve for the GFA. It can be seen that values of the pressure p and 
binodal and spinodal temperature of the system obtained within GFA 
are smaller than those obtained within MFA. 



4 Conclusions 

In this paper we have considered the M-component Ising model with 
site disorder. Several approximation for this model have been formulated 
here. Namely, we have obtained the expressions for the thermodynamical 
potential, free energy and pair correlation functions within TTA, GFA 
and MFA. 

The numerical calculations were performed for some simpler models: 
one-component Ising model, non-magnetic binary alloy and lattice gas. 
For example, we have chosen the nearest neighbours interaction. For the 
one-component Ising model we calculated the free energy as a function of 
the magnetization, temperature dependences of the magnetization m{T) 
and the two-tail A^^) (fluctuating parameter) within TTA, GFA, MFA. 
For non-magnetic binary alloy we obtained phase diagrams (binodal and 
spinodal curves) within GFA and MFA. For lattice gas the isoterms and 
coexistence curves are depicted. 

This work was supported in part by the International Soros Science 
Education Program (ISSEP) through grant No. PSU062015. 
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Figure 5: The p — V diagram for the lattice gas by various temperature 
(sohd curves - two-tail GFA, dashed curves - MFA; 1 - binodal, 2 - 
spinodal) 
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